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Анотація: Стаття розкриває стан сучасної наукової думки до питання розрахунку 
кількості пасажирів на авіаційному транспорті. Дослідження має оглядових характер, задачею 
якого є визначення Й оцінка методів при вирішенні задач транспорту. Для досягнення 
поставленої задачі використано наукові публікації, які уключено до наукометричної бази 
Scopus. Відповідні наукові підходи до визначення параметрів пасажиропотоку на авіаційному 
транспорті досліджено у розділі частини обрання та використання методів для розрахунку 
відповідних параметрів. Висвітлено авторську думку до питання обгрунтування вибору 
використаного авторами методу та оцінку адекватності запропонованих авторами методів. В 
роботі було використано метод системного аналізу при вирішенні поставленої задачі із 
проведення дослідження стану сучасної наукової думки про обрання методів розрахунку 
пасажиропотоку на авіаційному виді транспорту. 

Ключові слова: пасажиропотік, авіаційний транспорт, моделювання, пасажирські 
транспортні системи. 


1. Введення 


Питання планування параметрів транспортних систем, як й питання управління такими 
системами є актуальним питанням сучасної наки та практики. Відомо, що перешоджерелом 
розмірів й складу таких систем є пасажиропоток, яким визначаються такі параметри як: 
потужність автостанції, тип й кількість рухомого складу, інтервал руху, рисунок маршрутної 
мережі та інші. Відповідно, при проектуванні окремих елементів пасажирської транспортноїї 
системи мається необхідність розрахунку відповідного значення й характеристик 
пасажиропотоку у часі. Цим зумовлено наукові інтереси до визначення значень 
пасажиропотоку. У підтвердження актуальності визначеного питання виступає й значна 
кількість наукових статей, які присвячено таким дослідженням. 


2. Об'єкт та предмет дослідження 


Об'єкт дослідження - пасажирська авіаційна транспортна система. 
Предмет дослідження - методи визначення параметрів пасажиропотоку в часі. 
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3. Ціль та методи дослідження; 


Ціллю дослідження є визначення сучасних методів моделювання параметрів 
пасажиропотоку, опис таких методів й їхня оцінка. Для досягнення поставленої в роботі цілі 
було використано методи системного аналізу дослідженої літератури. 


4. Аналіз літератури. 


Зростаючий попит на послуги аеропортів цивільної авіації створив у світі актуальним 
питання визначення параметрів даного попиту. Накопичені бази даних параметрів 
експлуатації портів та пасажирських перевезень, демонструють значні відмінності між 
окремими аеропортами в одній й той самий час. Цим обумовлено потребу визначення 
параметрів експлуатації аеропортів та засобів транспорту авіаційної галузі у часі. Визначення 
таких параметрів на певний проектний час забезпечить процеси планування діяльності 
окремого порту у певному часовому інтервалі. Це є необхідним для створення стратегії роботи 
аеропортів із задовільним рівнем якості обслуговування пасажирів та організації 
експлуатаційних заходів засобів транспорту. Подібні питання й підходи до їхнього вирішення 
вже було висвітлено в роботах сучасників [1-4]. Авторами роботи [4] запропоновано підхід до 
визначення пропускної здатності аеропорту та коливань пасажиропотоку за допомогою 
введеного ними індексу авіапасажирів. 

Індекс авіапасажирів (АРІ): Встановлення значення X; як пасажиропотік аеропорту за 
одиницю періоду, АРІ для цього періоду визначається як Х, : 


ж Х =; X й 
X, — t min (1) 
X rix Е X min 


де Xmin Ta Xmax - Мінімальна та максимальна кількість авіапасажирів за одиницю часу 
відповідно, тоді як X й коливається від 0 до 1. 


Рівень індексу авіапасажирів (1.АРІ): Набір { Pis Pa, P3- р,} це набір послідовності АРІ з 
кількох одиниць часу. Після кластеризації створений кластер {N,} це сукупність об'єктів 


даних. Коли період АРІ становить 1 місяць, і місячний АРІ аеропорту складає Х,, рівень 
індексу авіапасажирів (1.АРІ) виводиться як: 


І, р, є (0,ї) 
2, р, є (i, j) 
N =43,p,€(j,k) (2) 


N, p, € (0,1) 


де Pi, рг, р: € одиницями часу, тоді як i, j, k, 0 є граничними значеннями кластера. 
Інформаційну ентропію можна використовувати для вимірювання ступеня невизначеності 
системи (або ступеня впорядкованості). Виводиться за такою формулою: 


H(X)= -Y P(x, log, Р(х) (3) 


ізі 


де P(xi) - ймовірність вибірки Xj, і п - кількість зразків. 
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Можна помітити, що чим менше ймовірність виникнення події, тим вищі значення 
інформаційної невизначеності та ентропії. Авторами запропоновано прийняти спільний 


розподіл ймовірностей випадкового вектору (X,Y) буде pij, потім двовимірна спільна ентропія 
вектора (X,Y) є: 


H(X. Y)=-$ У ру log p, (4) 


i=l j=l 


Вважаючи, що спільні розподіли ймовірностей X i У є pij та Pej, відповідно, умовну 
ентропію можна визначити як: 


H(X/Y)=- у У ров Ру 

i=l j=l Py (5) 
H(X/Y)=- \ У p, log Pi 

i=l j=l Pig 


Таким чином, можна виразити значення ентропії, для якого змінна X (або У) зменшується 
через появу змінної У (або Х). 


1(X;Y)=H(X)-H(X|Y) 
=H(Y)-H(Y|X) (6) 
= Н(Х)+Н(Ү)-Н(Х,Ү) 


Комбінуючи формули (3)- (6), повний вираз можна звести до такого вигляду: 


Pi 


I(X;Y) =) p; log, (7) 
i,j 


ig F gj 


Регресія опорного вектора є широко використовуваним методом прогнозування. Загальну 
модель лінійної регресії можна виразити так: 


ҒО) = м'х+Ь (8) 


де w - вектор нормалі вхідного вектора АРІ, i b— значення відхилення. 
Значення дорівнює нулю лише тоді, коли /(х) є точно таким же, як справжнє значення. 
Таким чином, концепцію можна виразити так: 


min Ww? ЄМО СЮ - уд.С»0 (9) 


i=l 


де С є константою регуляризації для виконання компромісного розрахунку на передній Ta 
задній панелі. 


0 if |z] «є 
С) = 10) 


|z] —é otherwise 
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У фактичних даних API певне значення може перевищувати звичайну тенденцію через 
зовнішні причини. Тому в разі серйозних відхилень від фактичне значення, слабкі змінні &, 


та Є вводяться як інтервали «пом'якшення», що зводить формулювання до наступного 


вигляду: 


тіп К н + cy (2 E) 
i=l 


т Ро) у, 5+8, 
Е у- (х) <е+6 


(11) 
i=1,2,...,m;č, é 20 


Використовуючи подвійний принцип і вводячи Лагранжів мульти-наконечники о; Ta а; , 


подвійну проблему SVR можна сформулювати як: 


max У у (а -а)-- 
пах 2.3 годе 

Ума; -aa чад x, (12) 
i=l j=l 


st.) (а; - а,) є 0,0 <a;,a,<C 


i=l 


Коли передбачене значення API потрапляє в є-м'який зона, а; та а, може бути 


ненульовим значенням. Нарешті, функція прогнозування регресії SVR може бути виражена 
як: 


(х) = Уа; задх'хчь (13) 


i=l 


b=y,+e-) (а, –а)хтх, (14) 


j=l 


Для даних часового ряду API з нелінійною тенденцією SVR може відобразити вибірку у 
просторі високої розмірності через функцію нелінійного відображення ф(х) , а потім замінити 


внутрішній векторний добуток простору високої розмірності ф(х,)Ф(х,) з функцією ядра 


K(xi,xj). Найбільш часто використовуваною функцією ядра є функція ядра з радіальною базою 
Гауса (КВЕ), яку можна виразити таким чином: 


K(x,,x,) = ех be-all 
із X;) = ехр( F ) (15) 


де гамма є параметром функції ядра Гаусового радіального базису 


(гамма = ) Ta o > 0 – пропускна здатність ядра Гауса. 


2 
С 
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Зрештою, функція регресії набуває такого вигляду: 


т 


f(x) = У (a; -а,)К(х, х) +b (16) 


Щоб проаналізувати результати прогнозування різних моделей, у цьому дослідженні 
використовували середню абсолютну відсоткову помилку (МАРЕ) 1 середньоквадратичну 
помилку (RMSE), які можна отримати за допомогою таких рівнянь: 


п 


2 


МАРЕ = =" x 100% (17) 


у, ЗУ; 


(18) 


де yii у, - це фактичні та прогнозовані значення. 


Запропонований авторами метод обчислення прогнозованих значень, в свою чергу, 
опирається на досягнення попередніх досліджень та може бути використаний при рішенні 
задач з організації перевезень пасажирів в частині станційного обслуговування та планування 
станційної діяльності. До переваг запропонованого методу можна віднести можливість 
планування навантажень на станційні зали очікування, місця стоянки засобів транспорту та 
їхній рух територією авіаційного порту. Недоліками запропонованого методу є 
короткострокове планування, неврахування коливань пасажиропотоку та відсутність в моделі 
інших видів транспорту. Відомо, що авіаційний транспорт, у певній мірі, може міняти 
розклади руху в режимі реального часу, що здебільш пов'язано із переносом рейсів до 
настання льотних погодних умов. Однак, такі ризики авторами не розглянуто й не введено до 
запропонованої моделі. 

В умовах зміни розкладів руху пасажир має можливість у користуванні не лише іншими 
авіаційними рейсами (альтернативними маршрутами), а й має можливість обрати інші види 
транспорту. Такі дії пасажирів обумовлені підсвідомим вибором між варіантом очікування 
настання сприятливих умов для здійснення авіарейсів та реалізацією потреби в переміщенні 
іншим видом транспорту із більшим часом їздки. При таких процесах вибору пасажир 
підсвідомо обирає, як для нього комфортніше й ефективніше здійснити їзду. З аналізу даного 
питання можна стверджувати, що у такому випадку питання комфорту та ефективності часто 
є похідною від часу, а вибір пасажира є підсвідомим у випадках наявності можливості у такому 
виборі. 

Авторами робот [5-6] розглянуто імовірність вибору пасажиром маршруту в межах одного 
виду транспорту. Моделюванням передбачено, що перевізні можливості маршрутів туди й 
назад майже однакові, про що свідчить про збалансованість руху авіапасажирів. 

У запропоновані авторами моделі введено обмеження, що вірогідність пересування одним 
маршрутом з міста А в місто В є однаковою (Р, , = Р, ,). Маршрут АВ приваблює пасажирів 
з міста А i міста В, у тому числі маршрут AB конкурує з тими маршрутами, які включають 


міста А i місто В. Вірогідність вибору маршруту (Рал-в) запропоновано розраховувати за 
залежністю 19: 
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C 
P =P = ae (19) 


A-B B-A 
ієМ№ д {+В ~i,A jeNg,j+A ~ Ј,В 


де Са,в - позначення конкурентоспроможності. 

Чим розвиненіша економіка двох міст, тим більшою є їхня привабливість до пасажирських 
кореспонденцій між ними. Авторами зазначено, що конкурентоспроможність маршруту 
пов'язана із рівнем розвитку конкретного міста A та В. Введено параметри Fa та Fg — 
привабливість міст A та В, відповідно. Запропоновано спосіб визначення Сдв - позначення 
конкурентоспроможності за допомогою рівняння 20: 


Cap =F, F (20) 


Привабливість окремого міста, Ha думку авторів, є похідною від кількості його населення 
та рівня економічного розвитку. Цим обумовлено припущення про можливість визначення 
привабливості міста використовуючи параметри кількості мешканців та рівня внутрішнього 
валового продукту (ВВП) в конкретному регіоні. 

Недоліком такого ствердження можна назвати неврахування соціального розвитку в 
регіоні, попереднім аналізом літератури визначено, що даний параметр впливає на параметри 
пасажирських кореспонденцій відносно окремого регіону. Визначено, що місця із розвинутим 
рівнем соціального стану суспільства мають більшу привабливість для утворення культурних 
та туристичних їздок до цих місць. 

Визначення пасажирських кореспонденцій між парою місць /(Та,Тв) запропоновано 
здійснювати за рівнянням 21: 


ІТ, Ту =a +bT y +CT.,;, (21) 


де Tmax ПОЗНачає пропускну здатність міста з більшою привабливістю; 

Tmin позначає пропускну здатність міста з меншою привабливістю. 

Авторами запропоновано використовувати метод найменших квадратів, щоб дізнатися 
значення а, b іс. 

У запропонований спосіб авторами запропоновано розрахунок пасажирських 
кореспонденцій між містами опираючись на технічні параметри аеропортів 1з застосуванням 
коефіцієнтів до таких параметрів. Однак, запропонований підхід не повною мірою відображає 
ствердження самих авторів про вплив стану розвитку регіону та кількість населення у певному 
пункті прибуття чи відправлення. Цим запропоновано вважати, що пасажирські 
кореспонденції залежать виключно від характеристик технічних можливостей самих 
аеропортів. Авторами не запропоновано методи визначення введених у залежність 2.21, 
їхнього калібрування та перевірки на відхилення. 

Запропонована авторами модель розрахунку імовірності вибору певного маршрут між 
пунктом А та В «Та-в» описана в залежності 22: 


<1 >= f ToTg) Pg (22) 


Запропонована модель враховує, що кількість пасажирів на маршруті не залежить від 
характеристик їздки на певному авіаційному рейсі. Не урахування таких певних чинників 
їздки, наприклад, час, комфортність та вартість такі розрахунки не є коректними, що 
визначено у попередньо описаних дослідженнях. 

Запропонована модель відображає процес вибору індивідуальних місць призначення: 
якщо місто А є привабливішим за інші міста, воно має більшу перевагу перед іншими містами. 
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Коли привабливості двох міст маршруту більші, конкурентоспроможність маршруту сильніше 
і пропускна здатність маршруту стає більшою. 

Одночасно, авторами запропоновано розраховувати пасажирські кореспонденції між 
містами і та j - Tij за запропонованою залежністю 23 із використанням гравітаційного 
моделювання: 


parha 0. (23) 
ое Ду 

де т; - населення міста, де розташований аеропорт і; 

пу- населення міста, де розташований аеропорт j; 

dij - відстань між містом, де розташований аеропорт і, та містом, де розташований 
аеропорт j; 

b=-0,13. 

Запропонована авторами модель 23 враховує вплив кількості мешканців у місцях і Taj та 
не враховує ані соціального ані економічного стану розвитку такої пари місць. Одночасно 
авторами обрано відстань перельоту фактором опору для реалізації таких пасажирських 
кореспонденцій, а комфортність, час та вартість даного переміщення не враховано, що 
неповною мірою розкриває фактори протидії кореспонденції. 

Можна стверджувати, що авторами в роботі не повною мірою враховано вплив на 
кількість пасажирів окремо розглянутого маршруту, як час, вартість комфорт, інтервал та інші. 

Автори робіт |7-17| також розглядали питання розрахунку параметрів пасажирських 
перевезень авіаційним транспортом. В роботі (10) авторами було викладено свій підхід до 
розрахунків певних параметрів із використанням теорії гравітації. 

Модель гравітації розроблена на основі відповідного закону, де вона обчислює силу 
тяжіння між двома об'єктами. Сила залежить від таких факторів тяжіння, як маса і відстань. 
Чим більше і ближче об'єкти, тим більше сила між ними. Запропонована гравітаційна модель 
для обчислення взаємодії або тяжіння між двома географічними місцями, такими як валовий 
продукт та населення. В економіці авіакомпаній гравітаційна модель, як показано в 24, 
використовується для прогнозування кількості пасажирів. Фактори привабливості впливають 
на попит на авіаперевезення. Приклади факторів впливу з місця розташування (A) - це ВВП i 
населення. Параметри, що стосуються маршруту (P), такі як відстань, час у дорозі або вартість, 
можна використовувати як фактори, що впливають на маршрут. Чим більше вплив впливають 
факторів, тим більше пасажирів подорожує між цими двома пунктами. 


P, =c] [AA p] [R (24) 
k=1 k=1 


де Pij— кількість пасажирів між пунктами і Ta j; 

Aike k” впливаючий фактор від місця розташування і; 

Rix є К" впливаючий фактор від маршруту між локаціями іта /; 

п і т відповідно кількість факторів, що впливають на розташування та маршрут; 

с є константою для всіх місць і маршрутів; 

Gk 1 ук є впливовими вагами місцеположення та факторів, що впливають на маршрут 
відповідно. 

Фактори впливу в 24 розділені, як і в 25, так що фактори впливу в обох кінцевих місцях 7 
та j мають різну вагу впливу, тобто ax i fx. 


P, = f Jas TART TR.” 9) 
k=1 k=1 k=1 
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У цій роботі 24 1 25 використовуються як перші дві моделі для оцінки потоку 
авіапасажирів. Запропонована Й модель 26: 


F= с] АХ П R” (26) 


Запропоновані авторами підходи не містять розкриття значень та смислу деяких 
складових в запропонованих моделях. 

Дослідниками [10-17] розкривалось питання вивчення поведінки пасажирів у питаннях 
повторного вибору певного способу пересування. Так, в роботі [17] описано ймовірність того, 
що авіапасажир вибере певну авіакомпанію Pro(t) й запропоновано залежність 26: 


чі 


Рго(т) = УРго(г |1) = Уе " (26) 


де t це поточна дата; 

1-10 використовується для опису періоду між останньою датою подорожі та поточною 
датою; 

Т - представляє період максимального інтервалу часу; 

ti— проміжок часу між найближчими подвійними датами подорожі; 

itait+l. 

Pro(t) - ймовірність того, що авіапасажир вибере цю авіакомпанію. 

Відповідно до запропонованої 26 залежності чим ближче між часом f і часом гарячого 
циклу, тим більша ймовірність вибору цієї авіакомпанії. Pro(t/;) використовується для 
позначення ймовірності вибору авіакомпанії в момент часу t під впливом гарячого циклу fi. 

Також, авторами запропоновано алгоритму розрахунку параметрів авіаперевезень 
пасажирів. Припущено, що Pm - матриця внутрішньої рушійної сили пасажира авіаційного 


транспорту, то відповідно можна описати шлях як PA. Запропоновано матрицю 
розрахунків у вигляді залежності 27, яка враховує таки параметри впливу внутрішньої 
рушійної сили пасажира, сили впливу інших пасажирів і подібної сили впливу авіаційного 
пасажира, які можуть бути розраховані за допомогою: 


R,,, = аР, +a,C,,,P,,, +0- –0,)8, Р. (21) 


пт nn nm 


де а, означає внутрішню рушійну силу повітряного пасажира; 
о, представляє тих авіапасажирів, які постраждали від інших пасажирів; 
1-а, - а, описується як пасажири, на яких постраждали подібні авіапасажири. 


Коефіцієнт відкликання р використовується для опису співвідношення, що прогноз 


авіакомпанії правильний. Для набору результатів прогнозування /, Ро; використовується для 
представлення потенційних наборів авіакомпаній, які можуть вибрати пасажири, і набору 
авіакомпаній, які фактично вибирають А; для авіапасажирів р; отримано з експериментальних 
даних, то коефіцієнт відкликання можна розрахувати за допомогою 28 та 29: 


пит(А, яф and A, eP, 
Г пит(А, # ф) 
num(A,#@ and А, €P, 
пит(А, + ф) 


(28) 


(29) 


International Science Journal of Engineering & Agriculture 2022; 1(3): 247-272 255 


де Prec - це рівень точності прогнозування системи; 
Ас представляє коефіцієнт точності системного прогнозування, який означає, чи 
подорожує авіапасажир чи ні, правильне прогнозування займає загальне прогнозування згідно 


з 30: 


А num(A,#@ and А єС and A, cE 30 

і пит(А, + ф) (30 

Якщо Fg - середнє гармонійне число коефіцієнта відкликання та швидкості точності, i 8 

це рівень важливості для коефіцієнта точності, що займає коефіцієнт відкликання, тоді 
співвідношення запропоновано авторами розглянути у вигляді 31: 


со 2 Вр 31 
Е авер 


Побудована авторами матриця результату вибору авіакомпаній Mim, де Mi 
використовується для представлення того, чи є авіапасажир р; вибирає авіакомпанію а), i 
значення дорівнює 1, поки воно відповідає дійсності, інакше значення встановлюється на 0). 
Оцінка сортування Rscorei визначається як 32: 


т 


К, = М, 2 К, аку (32) 


ў 
j=l 


I точність сортування 7) можна отримати шляхом запропонованим авторами у рівнянні 33: 


Е num(RS, #0) (33) 


Dis RS; 


Вплив кількості пасажирів Ha параметри перевезень є очевидним та впливає на доцільність 
функціонування всієї транспортної галузі. Авторами робіт [17-20] обчислювались питання 
доцільності функціонування аеропортів з підходу відповідності грошових потоків до 
економічної доцільності утримання такого елементу інфраструктури галузі транспорту. 

Авторами [20] обчислено беззбитковість пасажиропотоку для регіональних аеропортів. 
Беззбитковість, на думку авторів це кількість користувачів послуг авіа порту, яка забезпечує 
фінансові надходження відповідні для функціонування порту. Для моделювання відповідної 
кількості пасажирів - користувачів послугами порту, авторами запропоновано математичне 
моделювання. Висунута авторами модель для відповідного розрахунку 34 є лінійною: 


ТС, = ЕС, «УС, x Pax, + £, (34) 


де ТС; - загальна експлуатаційна вартість аеропорту; 

ЕС; - постійним компонентом експлуатаційних витрат; 

УС; - змінною складовою експлуатаційних витрат, що залежить від кількості 
обслуговуваних пасажирів Рах. 

Авторами роботи не визначено значення всіх складових запропонованої функції, що 
призводить до неможливості її повного розгляду. Визначено авторами, що рівняння 34 дає 
взаємозв'язок між експлуатаційними витратами та кількістю пасажирів у вигляді залежності 


35: 
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TR, = P x Pax, чи, (35) 


де ТК; - загальний операційний дохід, отриманий компанією аеропорт. 

Використання теорії графів та графоаналітичний аналіз сучасними науковцями було 
використано при моделюванні найкоротших відстаней між аеропортом та адміністративними 
центрами певного регіону. В роботах [21 - 24]. В роботі [24] авторами запропоновано модель 
36: 


L, = Min(£,)G=1,...m; j - 1,.. п) (36) 


де і є територією округу; 

j це аеропорт; 

Lij- геометрична відстань між і та J; 

k є найближчим аеропортом до i; 

а Lik це найкоротша евклідова відстань. 

Кожному округу приписується внутрішній район найближчого аеропорту, незалежно від 
того, чи обслуговується він аеропортом у межах його власних кордонів. Одним із показників 
авіатранспортного сполучення між парою міст з точки зору пасажиропотоків є абсолютна 
інтенсивність сполучення (Ti), яку авторами запропоновано розраховувати відповідно до 
рівняння 37: 


а) (37) 


де Jj 1 Її повітряні пасажиропотоки з міста і до міста j 1 з j до i, відповідно. 

Щодо авіаційного пасажирського транспорту, lij в основному має дорівнювати І. Аналіз 
домінантного потоку також було використано авторами для розв'язання певних задач з 
моделювання потоків пасажирів. Серед усіх можливих потоків з одного конкретного міста та 
всіх інших міст у мережі аналіз домінуючого потоку визначає лише найбільший потік із цього 
міста або до нього, як наведено у рівнянні 38. 


1 


ТТ. 
L, =max äl (=1,2,3,..п ken) (38) 
0, +D 


i і 


де Ті; 1 Tj це потоки повітряного руху з міста і до j та з j до і відповідно; 

Oi i Di загальна кількість авіапасажирів з міста та в місто і відповідно; 

Lik є «домінуючим коефіцієнтом потоку» міста і. 

Серед багатьох зв'язків між і та j (/=1, 2,...п), зв'язок ік є «пріоритетним зв'язком», що 
містить домінуючий потік. При цьому, науковці [21-55] виявляли критерії формування потоків 
пасажирами відповідно до корисності кожного з варіантів переміщення. Авторами [52] 
запропоновано дві моделі дискретного вибору, які припускають, що існує шкала переваг перед 
альтернативами, i людина вибере альтернативу з найбільшою корисністю. Модель 39 є 
найпростішою та широко використовуваною моделлю дискретного вибору для розуміння 
поведінки людей: 


Up =V tE; (іє, (39) 


де Vnj € спостережуваною корисністю альтернативи J; 
Enj - це неспостережувана корисність (або термін помилки) альтернативи j; 
Jn є сукупністю всіх альтернатив, які індивід п може вибрати. 
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Друга запропонована модель 40 може допомогти зрозуміти, що віддає перевагу 
респонденту серед нескінченних альтернатив. Модель 40 передбачає, що умови помилки 
(необслуговувана корисність) між альтернативами однаково та незалежно розподілені. 
Припущення означає, що перевага між двома варіантами не залежить від інших варіантів у 
наборі вибору. У моделі функція корисності індивіда п вибір альтернативи ) серед усіх 
альтернатив подається так: 


UE Vora. (jeJ,) (40) 


де Vnj є спостережуваною корисністю альтернативи j; 

Enj— це неспостережувана корисність (або термін помилки) альтернативи j; 
Jn є сукупністю всіх альтернатив, які індивід п може вибрати. 
Спостережувана корисність У; альтернативи j виражається таким чином: 


К; 
Vij = at XB (41) 
k=l 


де oj є константа сприяє альтернативі | спостережувана корисність; 

Xnjk Є атрибутами або пояснювальними змінними, які можуть вплинути на корисність 
альтернативи ); 

Вк є параметром пояснювальних змінних; 

К; відноситься до кількості пояснювальних змінних, пов'язаних з альтернативою j. 

Отже, в моделі 39 ймовірність вибору індивідом п альтернативи j серед усіх альтернатив 
Jn показано таким чином: 


- XPV) | (42) 
: х, ied, exp(V,,;) 


Припущення 40, використане в моделі 41, ймовірно, буде нереалістичним у ряді умов. У 
моделі 41 функція корисності індивіда п вибір альтернативи j під гніздом Р дається як: 


О. = H pV, 


і п Oa (45) 
де ць є параметром масштабу, оціненим за даними, значення негативно впливає на 
дисперсію терміну помилки; 
Упїь Є те саме, що спостерігається функція корисності в моделі; 
Lub Vnjlb є фактичною спостережуваною корисністю альтернативи j під гніздом b в моделі. 
Функція корисності гнізда b обчислюється на основі спостережуваної корисності всіх 
альтернатив усередині гнізда, що показано в 46: 


1 
У, =4 | —In У exp( HV ip (46) 


ДЬ Јер 


де Ap - параметр масштабу, пов'язаний з гніздом Р. 

Отже, співвідношення Ль/идь буде нормалізовано до //идь, який називається як параметр 
інклюзивного значення або параметр logsum гнізда b фактично оцінено в моделі 40. Отже, 
ймовірність того, що індивід вибере альтернативу j під гніздом b у моделі 46 було розраховано 
за наступною формулою 47: 


258 Dolia Olena: Analysis of modern scientific approaches to calculating the number of passengers on air transport 


exp(z,,V.,,) 
Ey Pop = ai Рон 
De jep SPU Vip) 


1 (47) 
exp| — In (5 Е ехр( H oV nip ) 


1 
Donen ®XP окр) 


ЛЬ 


де Рдь = умовна ймовірність того, що індивід п вибирає альтернативу j в гнізді b; 

Pnb це Ймовірність того, що індивід п вибрав гніздо b в межах набору гнізд В. 

Спостережувана функція корисності моделі дискретного вибору для цього дослідження 
запропонована у вигляді 48: 


У, = Constant + В.ТтауеїСозі,, + B,AccessTime , + B,JourneyTime, (48) 


+ B, Frequency , + B,SeatComfortMid , + B,SeatComfortHigh, + 


де TravelCost; це вартість проїзду (ціна квитка або вартість водіння) для використання 
альтернативного режиму /; 

АссеѕѕТіте; – час доступу до автовокзалу чи аеропорту; 

JourneyTime; це час у дорозі від пункту відправлення до пункту призначення (години); 

Етедиепсу; - кількість діючих автобусів або рейсів на тиждень; 

SeatComfortMid; 1 SeatComfortHigh; є фіктивними змінними, що представляють середній i 
високий рівень комфорту сидіння відповідно. 

Щоб краще зрозуміти конкуренцію між автомобілем, автобусом і повітрям подорожі в 
регіональній, авторами розраховано пряму та перехресну еластичність для вартості подорожі, 
часу поїздки та частоти обслуговування результатів моделі 47, використовуючи дані 
респондентів авіапасажирів та не-повітряних пасажирів відповідно. Пряма еластичність може 
представляти відсоткову зміну залежної змінної (наприклад, ймовірність вибору опціону), 
викликану зміною на один відсоток пояснювальної змінної (атрибуту), що представляє 
інтерес. Пряма функція пружності моделі 47 показана як: 


Ey, =11- P) + uy, -DU - Py BX (49) 


де Xik представляє пояснювальну змінну K альтернативи і; 

Р; це ймовірність того, що індивід п вибрав альтернативу і; 

Ридь- умовна ймовірність того, що індивід п вибір альтернативи і в гнізді b. 

Для моделі 46, ць замість цього буде зафіксовано значення 1. Крім того, перехресна 
еластичність може вказувати на зміну ймовірності у відсотковому співвідношенні конкретної 
альтернативи в сценарії вибору через граничну зміну зазначеної пояснювальної змінної іншої 
альтернативи. Таким чином, це важливий показник, який може представляти конкуренцію між 
альтернативами. Функція перехресної пружності відображається як: 


ЕЎ =-PP Xu (50) 
ЕЙ =A + (u p -DP plpa (5р) 


Рівняння (50) обчислює перехресну еластичність альтернативи ) щодо відсоткової зміни 
пояснювальної змінної k, за умови і та j в різних гніздах. Рівняння (51) відноситься до 
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розрахунку поперечної пружності за умови і та j з того самого гнізда. Функція логарифмічної 
правдоподібності є важливим індикатором для оцінки моделі, яка є ймовірністю, пов'язаною з 
даними спостереження за вибором. Функція логарифмічної правдоподібності оціненої моделі 
задається як: 


5 


x,Y)= УУУ, n P (|8) (52) 


n=l 5=1 j=l 


ШАВ 


де В параметр атрибута х оцінюється за моделлю; 

М кількість респондентів; 

5 це набір сценаріїв вибору; 

Ynsj дорівнює 1, якщо альтернатива j було обрано, інакше буде 0; 

Риз- оцінена ймовірність альтернативи j обраний фізичною особою n серед усіх 
альтернатив у наборі вибору 5. 

Авторами роботи [56] визначено, що вони вирішували питання прогнозування попиту на 
авіаційному пасажирському транспорті. В своїй роботі вони представили залежності 53 — 54: 


па Эъ) = Р, ооа (52) (53) 
Соу(х,,х,) = COV X an Xsan) (54) 

Y, =4Y, +0 +.. +a Y, +E, (55) 

Ү = u+ Y,- 4)+aY, = и) +: +a, (Үү, = ш) E, (56) 


де, ар- коефіцієнт авторегресії; 


Y, зе, - бе, - 0808 (57) 
Үү = а}, ої Ж a Yp tE, (58) 
OE OE T O Eq 

VY, =(1- BX, =Y, -Y,, (59) 
VY, =(1-B)Y, =(1-2B +В?) -2Y +Y, (60) 
Ф ,(B’)9(B)(1— В" - В')?7, = 8, (В')0 (В) є, (61) 

п LA 2 
RMSE = yo yi) (62) 

i=l n 

MAPE У З (63) 
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A 


n 3 
MAD =>. (64) 
i=l У; 

При цьому авторами не розкрито питання фізичного смислу наведених залежностей та 
значень приведених в залежностях складників. 

В роботах авторів [57-59] проведено моделювання обсягів пасажирських перевезень 13 
використанням Microsoft Excel i Matlab. Авторами запропоновано математичні моделі 
відповідного параметру для глобальної мережі аеропортів М (/=1,...,№ та поодиноких 
аеропортів T (f=1,...,7), загально система аеропортів приймає розмір МТ. Загальна 
математична модель пасажиропотоку між аеропортами запропонована авторами наведено у 


65: 


x, =Q, BT +l. 


1,1—1 


J5 ere зо T A ox; F (65) 


HAX + Nokis + доби a вва E O58 55 + Ui, 


де Xi в рівнянні (1) позначає k;X1 вектор внутрішніх змінних для аеропорту; 
ж ж . . . . 
x, k, x1 вектор зовнішніх змінних, специфічних для аеропорту; 


gi вектор | x 1 (тобто скаляр), що складається з однієї глобальної змінної - глобальної 
ціни на нафту; 

г позначає лінійний детермінований час тенденція; 

ии це специфічний для аеропорту процес ідіосинкратичної помилки, який, як 
передбачається, є послідовно некорельованим і слабо залежним з розміром ki х 1; 

ai, В, Ti, Ai, 1 0i вектори коефіцієнтів 1 матриці, які мають бути оцінено. 


Або із допомогою класичних одновимірних еталонних моделей розрахунку 
пасажиропотоку 66 (pi, di, gi) модель для аеропорту і читається наступним чином: 


ФО) pass, = а, + З (и, (66) 


де g(L)i 91.) позначають поліноми з відставанням порядку Pii qi. 

Проблеми сезонних (часових) коливань пасажиропотоків описано в роботі [60], у ній 
вводиться сезонна модель екстремального градієнтного підвищення (SD-XGBoost), щоб 
впоратися з проблемою прогнозування 1 в цій моделі також встановлено тимчасову кореляцію 
для часового ряду Y=y0,y1; "yn 9 - сезонний період, а сезонну різницю першого порядку 
визначено наступним чином 67: 


Vy, = У, - У, (67) 


де У єрізницевий оператор; 
5 l 
ye: Y, - зворотна різниця першого порядку з періодами 5 для зсуву для обчислення 
різниці. 
Різницеву послідовність з розміром ковзного кроку | i сезонним періодом 5 авторами 
викладено у 68: 


VOY == n У Na) (98) 


За допомогою набору із п зразків і т наведених ознак, 
D=({(x,,Vy,)}(Dl=n.x, cO”, V®y є 1 у, щоб вивчити набір функцій, які використовуються B 


моделі, авторами мінімізовано наступну регуляризовану мету виразом 69: 
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І (б) УКУ, Vy) + KD (69) 


де l є диференційованою опуклою функцією втрат, яка вимірює різницю між 
прогнозованим значенням сезонної різниці vy? У, іціль У? Y; fk B k-th. 

На відміну від звичайних методів оптимізації, кожна ітерація навчається шляхом 
додавання на основі попередніх ітерацій. vy У, дозволяє бути передбаченням і інстанція t 


ітерацію, 1 її можна передбачити на основі (7-1)-го члена, f; дозволяє звести до мінімуму 
відповідні помилки, що наведено у 70: 


LO -SIV VS, 0) +90) о 


i=l 


Завдяки розкладанню в ряд Тейлора другого порядку, рівняння 70 авторами 
запропонована записати як 71: 


LO = Sg, fa) + hSDN) 


nie 71 
where g, Sa V y D), (70 


h,=@ 


(5250-0) 


ТОУ VN) 


де gi i hi- градієнтна статистика першого та другого порядку щодо функції втрат, і постійні 
умови, які потрібно видалити, щоб отримати спрощену мету на кроці /. 
Рівняння 71 можна виразити у вигляді 72: 


го =) Ig о) + (x14 T+ Уа 
il 2 2 fa’ (72) 


-> ÈR o+ Хна o; |+yT 


ee з . 
Для нерухомої конструкції q(x), оптимальна вага @, j запропоновано обчислити 


рівнянням 73 та 74: 


обі 
зразками. (73) 
h+a 
ієї; l 
2 
б 1 T DIE 
Поу Ір нн (74) 
@- 2 од 


Рівняння 74 можна використовувати як функцію оцінки для вимірювання якості структури 
4. Конструкції g перераховуються за допомогою жадібного алгоритму. Ключовою проблемою 
навчання є знайти оптимальну схему сегментації на основі формули та прийняти точний 
алгоритм. Припущено, що 1/1 Jr є екземпляри наборів лівого та правого вузлів після розбиття 
розщеплення задається рівнянням 75: 
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Е a Ore: роз Узі? | РА 
К |У hA У. +2 Ў.т +4 


А . Й S) л . š 2 
Після 75 значення різниці ye У,у прогнозується через послідовність сезонних різниць 


5 z Я ia 
VOY, майбутнє прогнозоване значення авторами отримано за допомогою оберненої функції 
сезонної різниці, що містить спостережувані значення та значення різниці. Обернена функція 
сезонної різниці наведена у 76: 


1+1 


УФ =V ys EVO is TER EV Y 


i=S (76) 
=(y, у Yo) + Os =y) t+ Ya Yas) 
З рівняння 76, прогнозоване значення у;+; пасажиропотоку за час 1+ / є 77: 
Ў [Evy чуб +Y +y tet Ysa) = (Yas РУ OT y,) (77) 
i=S 


i=t+2-S 


AR ET 


Коли періоди зміщуються для обчислення різниці назад, тобто сезонний період 5-1, 


УУУ - різниця пасажиропотоку між /-м місяць і перший місяць, тобто 
Voy =y,- y, Toa у, зу, — Yọ: Рівняння 77 авторами зменшено до рівняння 78 
(5 (5 
умі ($v yı +V Sajt (78) 
ізі 


Для операції n виведення т моделей на першому етапі з наведеними зразками, 
п È ko зе or 
D = {(x,, у,)}, у Модель множинної лінійної регресії, як показано авторами в 79: 
I= 


Р(х) = ох + OX. +O, Xn +b (79) 


mim 


де Х; це вхідні дані; 
(@;b) - коефіцієнт регресії. 
Мета оптимізації моделі наведена авторами, як вираз 80: 


(Ф 5Б')- ага, тіпЕ,„, =arg min(y - X(@;b)) (80) 


(ab) 


де X — матриця з розміром пх(т- 1) 
При цьому значення в останньому стовпці завжди встановлюється як І, тобто: 
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хі Я. ЗОМ 1 
X= Ха Хо ОХ 1 (31) 
1 
Xn Xn2 7 Xam І Д 


В роботі Tsui W., Ozer Balli H., Gilbey A., Gow, Н. (2014) [61] авторами запропоновано 
математичну модель для розрахунку пасажиропотоку певного аеропорту із урахування 
сезонного коливання визначеного параметру. Модель для розрахунку пасажиропотоку для 
періоду часу – ARIMA (р, d, q), що не характеризується певними сезонними коливаннями часу 
авторами запропоновано у вигляді: 


Ф(р)У“Ү = а + (д), (81) 


Для розрахунку пасажиропотоку у певні проміжки часу, що має сезонні особливості — 
SARIMA (р, d, 4) запропоновано модель 82: 


Ф(р)®(Р)У“У?Ү, = а +9 (4)0(О)є, (82) 


де Ф(р) поліноміальний несезонний процес упорядкування р; 
w(P) позначає поліноміальний сезонний процес порядку Р; 
(4) позначає поліноміальний несезонний процес порядку 4; 
О(0) позначає поліноміальний сезонний процес замовлення Q; 


аур l : : TARE ; 
УУ - позначає рівень різниця для несезонних і сезонних процесів, відповідно, У, 


позначає залежну змінну, яку потрібно прогнозувати, є; позначає час помилки, 1 а позначає 
константу. 

Авторами, рівнянням 83 задано індикаторну змінну - 5, яка може мати постійний ефект 
або бути змінною в часі: 


_ П, якщо t>7 (під час і після втручання) (83) 
0, якщо г « Т(до втручання) 


або, Р, індикаторною змінною, яка може мати тимчасовий ефект або імпульсну функцію: 


p= 1, якщо? = Т(при втручанні) (84) 
0, якщо? + Т(до втручання) 


Коли втручання вставлене в модель SARIMA в рівнянні (82), модель втручання можна 
записати так, як показано в рівнянні 83: 


Ф(р)о(Р)У*У?Ү = а з (Н(ФО(О)є, + x, (85) 


де х; позначає функцію реагування (тобто постійну функцію або тимчасову функцію) або 
суму функцій реагування. 

Сезонні зміни пасажиропотоку та його коливання упродовж тижня розглядали автори [62- 
651. Запропоновано використання підходу із моделювання параметру кількості пасажирів при 
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використані методу регресії часових рядів. При використанні визначеного підходу до 
моделювання пасажиропотоку у певному аеропорті авторами отримано модель 86: 


М G G G 
узбіч У В 5, + ХУ. У, +У ФУ, а 519. Уа tN, (86) 
m=1 2-1 g=l g=l 


де 6 є параметром лінійного тренду; 

В сезонний параметр; 

Ути Є сезонною фіктивною змінною. 

Визначеним рівнянням передбачено, що якщо дані щомісячні, то M=12, квартальними — 
M=4 i так далі. 

У, € фіктивною змінною ефектів календарні варіації; 

У.г: € фіктивною змінною за місяць до виникнення ефекти календарних варіацій; 

Vg,t+1 Є фіктивною змінною через місяць після виникнення ефекти календарних варіацій. 

Якщо ефекти календарних варіацій зроблені щотижня, то С-4. Якщо ефекти календарних 
варіацій внесені в щоденні, то С-30 і так далі. 

Загальний вплив календарних змін можна визначити на основі часового ряду сюжету; 

М, це помилка білого шуму. 

Модель 86 запропоновано для включення значущих екзогенних факторів, які 
враховуються в моделі 87: 

M G G G Ө (B)@..(B* 
Y= B+ BS us EVV u E Vaa +L а tee ON, (87) 

m=1 g=l g=l ф, (В)Ф, (В") 


g=l 


Першими етапами побудови моделі 87 є створення регресійних моделей часового ряду з 
тенденційними, сезонними та календарними варіаціями з метою отримання похибки. 

Штучна нейронна мережа або відома як нейронна мережа, спочатку розроблена для 
імітації роботи людського мозку, складається з низки взаємопов'язаних простих елементів 
обробки, які називаються нейронами або вузлами попередньо продемонструвала себе у якості 
сучасного способу із моделювання поведінки пасажирів. Це призвело до можливості 
урахування поведінкових настроїв при певному моделюванні. Архітектура даного 
моделювання описує вихідне значення (ў) і входи (уг, уг2,...,угр) з математичною формулою: 


q р 
ў = Dwi pl) У иту +" | +b° (88) 
i=l 


j=l 


h , z . Е зі x 

де W; Є вага і-го нейронного вхідного шару до j прихований шар нейрона; 
b’ : . v v о М 
; Є Упередженістю j- го нейрону в прихований шар (/-/,2,...4); 


0 $ l а : 
W є Bara ј-го нейрона від прихованого шару до нейрона у вихідному шарі; 
b° є зміщенням нейрона у вихідному шарі; 
h В tes . .. tes 
f ; Є функція активації в прихованому шарі за допомогою сигмовидної функції, тобто f 


(x) = (1+exp(—x))"; 

f° є активація у вихідному шарі з лінійною функцією f (x) = x. 

Модель опорної векторної регресії долає проблему регресії, яка намагається отримати 
найкращу гіперплощину, використовуючи принцип мінімізації структурного ризику, шляхом 
поділу даних 1 мінімізації відстані між гіперплощиною та даними. Функція регресії методів 
авторами надано у вигляді 99: 
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Р(х) =w ф(х) +b (89) 


де w - ваговий вектор; 


ф(х) є функцією, яка відображає нелінійною х з вхідного простору у високорозмірні 
функції простору; 
Ь є упередженістю; 


w i b — коефіцієнти мінімізації функції ризику, описаної в наступному рівнянні 90: 


куо) =CLL Op fOD ні (90) 


0, if |y; - О) 
|у, = Р(х) 


«є 


where L,(y,, f(x) = 


-є, інакше 


де Le є є - нечутлива функція втрат; 


С i є є параметри, які були визначені для отримання оптимальних глобальних результатів. 
Концепція функції втрат полягає в мінімізації значення згідно 13 91: 


Rw éE) = тіпси + CDG 48) 91) 


з обмеженнями: у, - Ww ф(х) -Б<=+&, 
у @(x)-b<er€, 

wlx, )+b-y,<e+é,i 

EE >0, i=1,2,..,n. 


Оптимізація за такими обмеженнями може бути вирішена за допомогою використання 
методу Лагранжа 92: 


А я д 1 2 4 А 
Щоб, 2 E aa, Ba В) =z] р +E | 


п 92 
-У а wolx,)+b-y, нені | oF 


i=l 


«Уа | у нро) b+ +E |- BE + BS) 


i=l 


з використанням методу Каруша-Куна 92 можна представити у вигляді 93: 


(ааг) = Уа, а) -E$ (a, чаї) 
2 2 (93) 

Low А А 

= 72.3 а, зада; - а )УКОх,х;) 


i=l j=l 


Функція ядра K(xi,xj) може бути виражено як внутрішній продукт ф(х)! ф(х). Однією з 
функцій ядра, яка найчастіше використовується, є Гаусова радіальна базисна функція, яку 
можна визначити у даному випадку виразом 94: 
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2 
K(x,,x) =exp baal (94) 


2 
С 


де ос" є параметром ядра. 
Цим обумовлено вираз 95: 


код Уа, —a,)K(x,,x) +b (95) 


При прогнозуванні даних часового ряду вхідними параметрами є відставання даних 
спостережень х = Гуг-1, уг-2,..., Мі-рі. 

Автори [65] запропонували гібридну модель прогнозування, що складається з лінійної та 
нелінійної моделей, сформульованих у виразі 96: 


Y=L+N, +e, (96) 


де І, - лінійна складова 1 Ni складова гібридної моделі. 

В роботах |66-67) розглянуто попит на пасажирські та вантажні перевезення авіаційним 
транспортом у міжконтинентальному сполученні при різних умовах таких перевезень. 
Авторами ротом запропоновано нейронний метод математичного моделювання параметрів 
попиту. Авторами зазначено, що у даному випадку використовуються нейронні мережі з 
зворотним поширенням для покращення точності прогнозування попиту. Нейронна мережа 
складається з вхідного шару, вихідного тару і зазвичай одного або кількох прихованих шарів. 
Кожен з цих шарів містить вузли, 1 вони з'єднані з іншими вузлами на сусідніх шарах. Вихід 
із даного нейрона обчислюється шляхом застосування передатної функції до зваженого 
підсумовування його вхідних даних, щоб отримати вихід, функцією підсумовування є 97: 


I,= >; xX, (97) 


Функція активності представляє результат функції підсумовування й наведено в 98: 


net, = I, (98) 


Нелінійна передатна функція має сигмовидну форму згідно із 99: 


pot (99) 


net ; 


1+ехр ” 


Оскільки вихідний діапазон сигмоїдної функції знаходиться в межах від 0 до 1, 
відображення даних є більш придатним для використання із меншими відхиленнями. Оскільки 
не потрібно коригувати шаблон даних можна запропонувати наступну функцію: 


2) р р ур. (100) 


де Xmax 1 Xmin— це максимум і мінімум змінних; 
Dmax 0,8) i Din (-0,1) є зазначеними максимумом і мінімумом. 


International Science Journal of Engineering & Agriculture 2022; 1(3): 247-272 267 


Вузли у вхідному шарі представляють незалежні змінні проблеми, включаючи відповідні 
соціальні та економічні параметри. Прихований шар використовується для додавання 
внутрішнього представлення нелінійних даних. 


5. Методи дослідження; 


Для досягнення поставленої в роботі цілі було використано методи системного аналізу 
обраної літератури. 


6. Результати дослідження 


В результаті проведеної роботи було проаналізовано шістдесят сім літературних джерел — 
періодичних наукових видань включених до наукометричної бази даних Scopus. 3a 
результатами аналізу дослідженої літератури можна зробити висновок, що з часом питання 
визначення параметрів пасажиропотоків на авіаційному транспорті набирало актуальності й 
висвітлювалось частіше рік від року. Науковці використовували для рішення задач такі 
методи, як математичне та комп'ютерне моделювання. Серед методів математичного 
моделювання можна відокремити наступні: теорія графів, нейронні мережі, гравітаційний 
підхід та ймовірнісний. 


7. Перспективи розвитку подальших досліджень сучасних наукових підходів до 
розв'язання питання із моделювання пасажиропотоків 


Відповідно до проведеного дослідження сучасних наукових підходів до розв'язання 
питання із моделювання пасажиропотоків можна припустити можливість обрання методів 
комп'ютерного моделювання пасажирської мережі із її подальшим обчисленням при 
використанні теорій графів та тяжіння. 

Відповідне моделювання має забезпечити можливість внесення в мережу додаткових 
вузлів або ланок мережі для моделювання параметрів перевезень при внесенні змін в мережу. 
Використання теорії графів призведе до обчислення матриць відстаней та максимальних 
потоків по ланках мережі, а теорія гравітації розподілить значення параметру прибуття чи 
виїзду відповідно по ланках мережі. 


8. Висновки 


Відповідно до результатів можна зробити такі висновки: 
1. Задача моделювання пасажиропотоку є актуальною задачею сучасної науки й 
практики; 
2. Математичне й комп'ютерне моделювання є сучасними методами розрахунку 
параметрів пасажиропотоків. 
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Abstract: The article reveals the state of modern scientific opinion on the issue of calculating 
the number of passengers on air transport. The research is of a review nature, the task of which is to 
define and evaluate methods for solving transport problems. To achieve the task, scientific 
publications included in the Scopus scientometrics database were used. Relevant scientific 
approaches to determining the parameters of passenger flow on air transport are studied in the section 
of the selection and use of methods for calculating the relevant parameters. The author's opinion on 
the issue of justification of the choice of the method used by the authors and the assessment of the 
adequacy of the methods proposed by the authors are highlighted. In the work, the method of system 
analysis was used in solving the task of conducting a study of the state of modern scientific opinion 
on the selection of methods for calculating passenger traffic on the aviation mode of transport.. 
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